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1 Introduction 

Higher-dimensional black holes and black strings appear as classical solutions of gravity 
theories in more than four dimensions. One central issue concerning their behaviour is the 
stability issue. Recently, the stability of higher- dimensional black holes and in particular of 
black strings, is under intense investigation. The stability of higher-dimensional black holes 
has been addressed adequately in [1]. Also, the understanding of the stability of black rings 
and in general of higher-dimensional black objects has also been improved considerably (for 
a review see [2,3]). 

In the case of black holes that are solutions of gravity theories with high order curvature 
terms [4-9], the stability issue is an intriguing one. The Lovelock theory [10] is a non- 
trivial extension of General Relativity which apart from the Einstein-Hilbert term also 
includes higher order curvature terms. In fact, the Lovelock theory contains terms only up 
to the second order derivatives in the equations of motion which makes the theory more 
tractable. In spite of that, the black hole solutions of Lovelock theory could be found 
mostly because they are highly symmetric (for a review see [11], [12]). The study of their 
stability requires the application of linear perturbation theory, which however, confronts 
with the high complexity of the Lovelock equations, since the perturbative terms break the 
simplifying symmetries of the background metric. 

In the simplest case of a Gauss-Bonnet theory which is a second order Lovelock theory, 
the classical stability of black hole solutions has been studied. The stability analysis under 
scalar, vector and tensor perturbations has been performed [13-15]. It was found that 
there exists a scalar mode instability in five dimensions, a tensor mode instability in six 
dimensions, and no instability in other dimensions. Recently the master equations of 
Lovelock black holes of scalar, vector and tensor perturbations were derived in [16-18] 
and it was shown that small Lovelock black holes are unstable in the asymptotically flat 
cases [19]. 

Higher order curvature terms are known to appear in string theory, generically introduc- 
ing higher derivatives in the metric [20]. This will inevitably lead to perturbative ghosts. 
However, the second order curvature terms are known to take the form of the Gauss-Bonnet 
combination [21]. This corresponds to the non-trivial second order Lovelock theory which 
has no higher derivatives in the effective string action, so we expect no ghosts to appear 
in second order. However, it was found that spherically symmetric vacuum solutions in 
the Gauss-Bonnet theory suffer from ghost-like instabilities and it was conjectured that 
these instabihties persist in all vacuum Lovelock solutions [22]. In fact it was shown that 
there is a limit, the Chern-Simons limit, in which the theory becomes strongly coupled and 
linear perturbation theory breaks down. In this limit it was also shown, that if there is a 
fine-tuning between the parameters, the two branches of vacuum solutions coincide with 
the Chern-Simons black hole solution which has maximum symmetry. Scalar perturbations 
of the Chern-Simons black holes were studied in [23]. 

Higher order curvature terms also appear in braneworlds and especially in codimension- 

2 braneworld models which were mainly introduced because there is a mechanism of self- 
tuning of the effective cosmological constant to zero [24]. However, soon it was realized 
[25] that one can only find nonsingular solutions if the brane energy momentum tensor is 
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proportional to its induced metric. To reproduce the effective four- dimensional Einstein 
equation on the brane one has to modify the gravitational action by the inclusion of a 
Gauss-Bonnet term [26,27] (for a review see [28]) or to regularize the conical singularities 
(see [29] and references therein). 

Recently black hole solutions of codimension-2 braneworlds have been found. A six- 
dimensional black hole localized on a 3-brane of codimension-2 was proposed in [30]. How- 
ever, it is not clear how to realize these solutions in the thin brane limit where high order 
curvature terms are needed to accommodate matter on the brane. Generalizations to in- 
clude rotations were presented in [31] and perturbative analysis of these black holes were 
carried out in [32-34]. 

In the case that there is a Gauss-Bonnet term in the bulk, black hole solutions were 
studied in a five-dimensional codimension-2 braneworld model [35] . Two classes of solutions 
were found. The first class consists of the familiar BTZ black hole [36] which solves the 
junction conditions on a conical 2-brane in vacuum. These solutions in the bulk are BTZ 
string-hke objects with regular horizons and no pathologies. The second class of solutions 
consists of BTZ black holes with short distance corrections. These solutions correspond to 
a BTZ black hole conformally dressed with a scalar field [37,38] and they have black string- 
like extensions into the bulk. Generalizations to include angular momentum and charge 
were presented in [39] . Also four-dimensional Schwarzschild- AdS black hole solutions on the 
brane were found, which in the six-dimensional spacetime look like black string-like objects 
with regular horizons [40] . The warping to extra dimensions depends on the Gauss-Bonnet 
coupling which is fine-tuned to the six-dimensional cosmo logical constant. The presence 
of the Gauss-Bonnet term in codimension-2 braneworlds has important consequences. Its 
projection on the brane gives a consistency relation [27] that dictates the form of the 
solutions. It allows black string solutions in five dimensions, and in six dimensions it 
specifies the kind of matter which is needed in the bulk in order to support a black hole 
solution on the brane. We note here that black string solutions with a Gauss-Bonnet term 
in the bulk are not possible in codimension-1 braneworlds [41]. 

The stability analysis of these static black hole solutions of codimension-2 braneworlds 
is interesting. They solve N-dimensional Einstein field equations with the presence of a 
Gauss-Bonnet term in the action. The symmetries of the solutions are /A"^ x /C""*^ ^, where 

is a maximally symmetric space while the space K, is axially symmetric. Therefore, any 
stability analysis of these static solutions, unlike the Lovelock black hole solutions which 
in general are spherically symmetric, has also to confront with the particular symmetries 
of the solutions. 

In this work we address the problem of stability of codimension-2 black strings with a 
Gauss-Bonnet term in the bulk and more generally of gravity theories with high order cur- 
vature terms that have time-independent solutions with symmetries other than spherical.^ 
In codimension-1 braneworlds the Schwarzschild metric on the brane was considered and 
its black string extension into the bulk [44] was studied. It was found that this string is 
unstable to classical linear perturbations [45] . A lower dimensional version of a black hole 

^ Black hole solutions with a Gauss-Bonnet term and with these kind of symmetries were discussed 
in [42] 

^The stability of uniform black string solutions in the Gauss-Bonnet theory was studied in [43]. 
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living on a (2+l)-dimensional brancworld was considered in [46]. A BTZ black hole on the 
brane was found which can be extended as a BTZ black string in a four-dimensional AdS 
bulk. Their thermodynamical stability analysis showed that the black string remains a sta- 
ble configuration when its transverse size is comparable to the four-dimensional AdS radius, 
being destabihzed by the Gregory- Laflamme instabihty [45] above that scale, breaking up 
to a BTZ black hole on a 2-brane. The stabihty of BTZ black string was also discussed 
in [47]. 

One way to study linear stability of higher dimensional objects with curvature correc- 
tions is to find the eigenvalues of the Lichnerowicz operator of a given perturbation [48] . 
This method has the advantage of being formulated in a gauge invariant way allowing the 
study of metrics with various symmetries. We will derive the Lichnerowicz equation in 
the presence of the Gauss-Bonnet term. We will show that a simple application of the 
modified Lichnerowicz equation to the case of D=6 spherically symmetric black hole so- 
lutions of Gauss-Bonnet theory can easily reproduce the known results of [13] for tensor 
perturbations for these black holes. 

We will subsequently use the modified Lichnerowicz equation to the study of linear 
perturbations of five- dimensional black string solution of codimension-2 braneworld model 
of [35], away from the Chern-Simons limit. We will show that the knowledge of the Lich- 
nerowicz equation can provide very important information on the stability analysis of a 
complex system of coupled differential equations that has to be solved. In the case of 
codimension-2 geometries, the black string can propagate in two transverse extra dimen- 
sions so intuitively, one expects that metric perturbations to have more severe stability 
problems than the conventional five- dimensional black strings propagating in one extra 
transverse dimension. 

Our analysis shows that for tensor perturbations the modified Lichnerowicz equation, 

due to the symmetries of the codimension-2 black string solution, can not give any informa- 
tion on the perturbed system, indicating that either there is no tensor propagating modes 
or we are in a strong coupling regime. For the vector perturbations we find a degeneracy 
of the modes on the bulk space having no dependence on the brane coordinates indicating 
an instability on the modes generated on the brane and propagating into the bulk. To cal- 
culate the scalar perturbations, we will apply the derived modified Lichnerowicz equation 
to the scalar part of the metric perturbations. We will show that, due to the fact that the 
metric of the black string solutions can be brought to a factorizable form, the results are 
the same as having studied the propagation of a scalar field in the background metric of 
the black string, by solving the Klein-Gordon equation. We find stability for the scalar 
modes by both methods. 

Another important information that the modified Lichnerowicz equation can give us is 
the behaviour of the theory in the Chern-Simons limit. As we have already discussed, in the 
Gauss-Bonnet theory there are two branches of spherically symmetric vacuum solutions: 
The Schwarzschild-AdS branch (known as the Einstein limit) and the string branch (known 
as the Gauss-Bonnet limit) [22]. Both branches coincide at the Chern-Simons limit where 
the theory becomes strongly coupled. As we will discuss in the following, in our case the 
Chern-Simons limit manifest itself as a prefactor in front of the modified Lichnerowicz 
equation. The limit where this prefactor goes to zero is an indication of a strong coupling 
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problem, signaling that linear perturbation theory breaks down. 

The paper is organized as follows. In section 2 we present the general formalism for linear 
perturbations and derive the Lichnerowicz equation. We further generalize this formalism 
by including the Gauss-Bonnet term and derive the modified Lichnerowicz equation where 
a source term appears, due to the presence of the Gauss-Bonnet term. In section 3 we 
apply these results to spherically symmetric solutions of the Gauss-Bonnet theory. In 
section 4 we calculate the scalar perturbations of the five-dimensional black hole solution 
of the codimension-2 braneworld model and we discuss its vector and tensor perturbations. 
Finally, in section 5 we conclude. 

2 The Lichnerowicz Equation 

In this section we will review the general formalism of linear metric perturbations and we 
will derive the D-dimensional Lichnerowicz equation [48] and the D-dimensional modified 
Lichnerowicz equation in the presence of a Gauss-Bonnet term. 

2.1 General formahsm for hnear metric perturbations 

Consider a D-dimensional Einstein manifold {B,g). The stability analysis can be reduced 
in finding a solution of an ordinary differential equation of a Schrodinger type, i.e., the 
Lichnerowicz equation. This amounts to determine the spectrum of the Lichnerowicz oper- 
ator on transverse traceless symmetric tensor fields of the manifold B. We will first expose 
the general formalism and then we will apply it to a manifold with a Gauss-Bonnet term. 
We start from the definition for the Riemann Tensor as in [49], 



and consider the following linear perturbation of a £)-dimensional metric background n 



where the Latin capital indices M, N take values in the D-dimensional space, and all 
unperturbed quantities are written as X. We also decompose the Ricci tensor and a D- 
dimensional energy momentum tensor Tmn respectively as 



(2.1) 



(2.2) 




Rmn + £ SRmn, 
Tmn + £ STmn- 



(2.3) 
(2.4) 



The Einstein's field equations are 



Rmn — 7: Rqmn — Tmn, 



(2.5) 



and using the trace they can be written as 



2 




(2.6) 



Rmn — Tmn — 7: 



2 D-2 
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The zeroth order Einstein equations (or background solution) are 

1 2 _ 

Rmn — Tmn — 2 £) _2 ' 

while the first order Einstein equations (or perturbed equations) are 

SRmn — ^Smn, (2-8) 

where 

SRmn = -2 (□ W + VatVm h - V^Vm Hkn - V^Vjv W) , (2.9) 
with OHmn — g^^^A^shMN, and 

5Smn = STmn — 2 £) _2 i^MN g 5Tkl + Qmn h Tps — Hmn ) . (2-10) 

It is easy to check that the Einstein equations (2.5) are satisfied for qmn — 9mn + £ huN 
and Tmn = Tmn + £^Tmn- To simplify the equations we choose the de Bonder gauge, 
namely, the traceless and the transverse gauge conditions, respectively given by 

g'^'^hMN^O, (2.11) 

V^/lMAr = 0. (2.12) 
Then, in this gauge the first order Ricci tensor (2.9) can be written in the following form 

6 Rmn = —- (n Hmn — 2 Rknml h^^ — Rm ^kn — R-n hMx) ■ (2-13) 

In vacuum {Tmn — 0), the Einstein equation reduces to Rmn — 0, whose zeroth order 
part, Rmn = 0, can be used in the first order part (2.13) to obtain the Lichnerowicz 
equation for linear perturbations in vacuum 

nhMN-2 Rkmnl /i"^^ = 0. (2. 14) 

When a cosmological constant is present, le., {Tmn = —^dQmn)-, the zeroth and first 
order parts of the Einstein equations (2.6) are respectively 

Rmn — —-^ogMN, (2.15) 
SRmn — —-^d hMN, (2-16) 

where Ai^ = §z|Ad. Using (2.15) in (2.13), the first order part (2.16) gives again the 
Lichnerowicz equation (2.14). 
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2.2 The modified Lichnerowicz equation 

We want to establish the equivalent to the Lichnerowicz equation in the case where a Gauss- 
Bonnet term is included in the manifold B. In five or in six dimensions the Gauss-Bonnet 
density term in the action is given by 

aCoB ^a{R^- ARkl i?"^^ + RklpqR''''''^) , (2.17) 

where «(> 0) is the Gauss-Bonnet coupling constant. Its variation includes the Gauss- 
Bonnet term aHuN in the field equations with Hmn being the Gauss-Bonnet tensor, 

Gmn — ocHmn — —^DgMN- (2-18) 
For convenience, we split this tensor in five terms HMNi{i — 1, 5) which are 

Hmni — -^Qmn ^GB, (2.19) 

Hmnq — —'^RRmni (2.20) 

Hmnz — +4^Rmk (2-21) 

Hmn4 — +^Rmln Rr^ (2.22) 

Hmns — —'^RmklpRn^^- (2.23) 

As for the linear perturbations, we decompose the Gauss-Bonnet tensor as 

5 

Hmn — Hmn + £ SHmn — ^ {HMNi + £ ^Hmn^ , (2.24) 

i=l 

and after some lengthy calculations, where we used the gauge conditions (2.11) and (2.12), 
we find that the first order Gauss-Bonnet tensor decompositions are 

SHmni = -^j^hMN^^GB — Qmn Rkl R + '^Qmn R^^ ^hKL — '^Qmn R^ 

-guN R'"'-'''^ hKR R\pQ + 2-gMN i?'"'^'^'' VqVl Hkp, (2.25) 

SHmn2 — '^h^^ Rrl Rmn + RO Hmn — 2^ R^mn^ ^lk 

— RR\ hMK — R R^^M hNKi (2.26) 

SHmNs = '^R\ { — ^ ^MK + '2R^mk'^ ^'LQ + R\ hML) 

+2R\, (-□ hKN + 2R\k'^ hLQ + R^K hiN) , (2.27) 

^HmN4 — '^R^^ {~R^KLN ^MP + M hxN — ^ K hMN — ^ N^ M hxL + ^ N^ K hML) 

+2Rkmln i-h'^'' R^'p - □ /i^^ - 2i?^^^« hQp) , (2.28) 
SHmns — -\-2R^^^ {V pV xhuL — '^ p^ uhKL) + 2Rj^^^ {V pV xhiN — '^ p^ Nhxh) 

—Rmklp R9^^^ Hqm + 4:Rj^^^ Rnkqp h\ — Rj^^^ Rqklp h^M- (2.29) 

The Einstein's field equations (2.18) can also be written as 

Rmn = a Hmn - Qmn - ^ ^ ^ Qmn {aH - D Ad) , (2.30) 
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where the zeroth order equations are 

Rmn = tt Hmn - Azj gMN - {olH - D Ad) . (2.31) 

Using (2.31) together with (2.13) the first order field equations can be written as 

□ Hmn — '^Rkmnl = « ^mjv, (2.32) 

where 

2 - KL 

^MN = —'2SHmn + — - Qmn 9 SHkl — — - Qmn Hkl h 

+H^hKN + H§ Hmk- (2.33) 

When we calculated the variation of the Gauss-Bonnet term, the harmonic gauge was 
used. However, there is still some gauge freedom which can be further gauged away. 
The second term can be written as 



q'^'-SHkl ^i2-D) h^'^'KMNR + (8 - 2D) R^'^'u Hmn + (4D - 10) R^'"^ R^^^^shpN 
-{D-2) R^^^'R^'^^shAP + ^'^''R'phLT + {2D - 8) R^^'^'V mV bHas 

(2.34) 

while the third term becomes 

Hkl = —R Rkl + '^Rkp Rl + '^^'^kpl ~ 2RkmnpRl^^^ 

(2.35) 

Equation (2.32) is the modified Lichnerowicz equation, the difference with (2.14) is the 
appearance of the source term oBmn due to the Gauss-Bonnet term in the action. 

Another way of arriving to the same result is to vary directly the equation (2.18). 
However, the modified Lichnerowicz equation (2.32) is useful because it gives directly the 
perturbation equations once the form of the perturbation is known. 



3 Perturbation analysis of spherically symmetric Gauss-Bonnet 
black hole solutions 

As an application of the formalism developed in the previous section we will examine 
the perturbations of spherically symmetric solutions of Einstein-Gauss-Bonnet equations 

Gmn — ocHmn — —^D guN, (3.1) 
where Gmn is the Einstein tensor and Hmn is the Gauss- Bonnet tensor 

Hm = ij^QM ~ ^^klR^^ + RklpqR^^^^) — 2RRm (3.2) 

,A p TfNK I ^ p N rfKP r, p TjNKLP 

+4/£mK-K ^^KMP ^ —^J^MKLPJ^ 
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Namely, we will investigate the following equation 

5G^ - aSH^ = -AdS9a- 
In the transverse traceless gauge the variation of the Einstein tensor gives 

1 



(3.3) 



SGab 
and 



2 ^^^^^ ^ B^LK — R A^BL — R B^Al) — -hAsR 



SGa = q^^SGak - H^^Gak- 



QABh^'-RKL (3.4) 



(3.5) 



Furthermore, 

5Hf = g''''5HAK - h^^'HAK, (3.6) 

where the variations 5Hak are given by the equations (2.25)-(2.29) which already incorpo- 
rate the transverse and traceless gauge. Additionally, it is easy to see that 



5g^ =g''''hAK-h'''^gAK = 0. 



BK , 



(3.7) 



We will consider spherically symmetric black hole solutions of the Einstein-Gauss- 
Bonnet theory in D = 6 dimensions. The procedure is the same in the case of five di- 
mensions. The resulting equations are 



SG'a"" + a' SH'a"" = 0, 



(3.8) 



where SH'a^ are given by (2.25)-(2.29) and they coincide with the equations derived in [13] 
for tensor perturbations with the identifications —\Hab — H'ab and a — 2a' . 

Following [13] these equations can be solved transforming them to a Schrodinger-like 
equation. Consider the following spherically symmetric metric 



1 



-f{r)dt^ + -FT-rdr'^ + [de"^ + sin^ 6 {dip^ + sin^ (p {dx^ + sin^ X#^))] 



which satisfies equations (3.1). Perturbations of the above metric read [13] 

5'a/3 gap + hal3: 

where 



(3.9) 



(3.10) 



hap = (r, t) 








/ 

hee he 



(3.11) 



\ 



hex 

hipx h^ 
h^e hx<p h^x 

y h^Q h^ip h^x hipip ) 

For every mode in the perturbation we have hij — {9, (p, x, V') with i, j — 9, (p, x, V'- 
thermore, the transverse and traceless choice for the gauge alongside with the symmetries 
of the metric implies that the restriction of has to the sphere is transverse, traceless, and 
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so it can be expanded using a basis of eigentensors of the Laplacian on the sphere [50] . So 
with i,j running on the sphere we have 

hij^r^(f){r,t)hij{x), (3.12) 

where 

V^Vjiij^lhj, V% = 0, g'%j^O. (3.13) 

The bar refers to metric and tensors on the -S"^. 

Substituting (3.11) into (3.5) and using the expansion of hij according to the basis of 
eigentensors of the Laplacian on the sphere we find equation (11) of [13]. But somehow this 
is expected. In addition, we find the same equations as in [13] for the Gauss-Bonnet com- 
bination. Substituting again (3.11) to (3.6) and using both the expansion of h^j according 
to the basis of eigentensors of the Laplacian on the sphere, alongside with the transverse 
and traceless gauge, we end up after a lot of algebra to the equation (12) of [13]. 

The crucial point in the above analysis comes from the symmetries of spacetime. Here 
spherical symmetry allows for the decomposition of hij according to the basis of eigentensors 
of the Laplacian on the sphere, which at the end is the key factor for obtaining the final 
master equation of the perturbation analysis. Such a simplification is more difficult and 
some times not possible in the cases where the spacetime allows for symmetries other than 
spherical as we will discuss in the following sections. 



4 Perturbations of codimension-2 black strings 

We consider the following gravitational action in five dimensions with a Gauss-Bonnet 
term in the bulk 

(4.1) 

We are looking for solutions of the form 

dsj = g^y{x, p)dx^dx'' + a^{x, p)dp^ + L'^{x, p)de'^ . (4.2) 

The topology of the two-dimensional space (p, 9) can be represented by a cone of deficit 
angle (5. Regularization of this space dictates the introduction of a brane located at the tip 
of the cone. As we will discuss, the presence of the brane has important consequences in 
the perturbative analysis of these spaces. Solutions of the form (4.2) have been obtained 
in [35] . In this section we will derive the general formalism of metric perturbations of these 
solutions and we will discuss the scalar, vector, and tensor perturbations. 

4.1 General formalism of metric perturbations 

Gonsider the metric 



2 



-n'^{r)df + 4^ + r'd(j)'' 



+ dp^ + b\p)d9^ , (4.3) 
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which is of the form (4.2). In what follows we will consider /(p) = cosh(p/2y^), b{p) = 
2 (3 y/a sinh {p/2y/a), and n^{r) = —M + r^//^, which is a solution of the action (4.1) with 
the inclusion of a brane boundary term [35] . 

In order to investigate the stability of the solution we need to adopt an Ansatz for our 
perturbations. It is possible to see that the above metric has 3 Killing vectors. We can 
adopt a Fourier expansion of the perturbation with respect to the coordinates t, (f), 9 and 
keep all the components of the perturbation as functions of r, p . We are going to consider 
perturbations which retain the angular symmetry on the brane as well as on the transverse 
space. Keeping the axial symmetry both on the brane and in the transverse space means 
that everything must be Lie derived by 9^ and do to zero [32]. In this way, we will work 
on the s — wave approximation of our system for the two angular coordinates. 

As a result of the symmetries of the spacetime we can split the perturbation into a 
purely two-dimensional transverse piece, a mixed transverse and three-dimensional piece 
and a purely three-dimensional piece. This can be represented by 



(4.4) 



where (//, u — t,r, 0) and {i, j — p,9). In the Kaluza-Klein spirit, these perturbations can be 
interpreted as scalar, vector, and tensor, respectively, with respect to the three-dimensional 
spacetime. We will consider the following Ansatz for our perturbation 

u -^nt f (r, p) hf,i (r, p) \ 

The above ansatz contains the maximum information concerning the perturbation of 
our system, while at the same time is consistent with the aforementioned arguments. 

Substituting (4.5) into the Lichnerowicz equation (2.32), we can see in a straight forward 
manner that, the h^0{r\ p) mode decouples from the rest of the modes straight away, while 
the ht(i, (r, p) , /i^^ {f-i p) ^-nd hp^ (r, p) modes are coupled together and the same happens 
also for hte ij-, p) , hre ij', p) and h^e {r, p) modes. We could set them to zero in order to 
examine the rest of the system, but for the moment we will keep them, since in any case 
these modes, for the s-wave approximation that we are considering, do not interact with the 
other modes. As a first step we will examine the behaviour of the scalar modes hee {r,p), 
hpp{r,p) and hpe{r,p). 



4.2 Scalar perturbations-solving the Lichnerowicz equation 

We consider first the scalar perturbations. Our aim is to solve the modified Lichnerowicz 
equation (2.32) in the background metric (4.3) using the perturbation Ansatz (4.5). Sub- 
stituting our Ansatz for the perturbation we can see that all the perturbation equations 
are governed by a prefactor, which has the form {P — 4a). We see that when we are close 
to l'^ Aa, we face a strong coupling problem. We will refer to this limit later on. 

Using (2.34) and (2.35) we can calculate the two equations for the scalar modes. These 
equations couple the scalar modes of the perturbation with the tensor and the vector modes. 
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However, they can be decoupled. Here we will summarize the results while in the appendix 
we present the technical details. The two scalar modes are given by 



h, 



pp 



100 



haa — 



y^sinh (^j 

P 



2 ^ \foL sinh 



u (r, p) . 



(4.6) 



(4.7) 



Separating the variables choosing u (r, p) — f (r) y (p), the functions / (r) and y (p) are 
given by the differential equations 



2v^ 



3 cosh 



P \ dy{p) 
dp 



+ ^foL sinh 



P 



■9^ V (p) 
d p^ 



with m being the separation constant 



+ my (p) = 0, 

(4.8) 



m — 



f{r) + 



\ ^ or ^ or^ 



PM-r^J ' ^ ' Pr 
The equation (4.9) can be written in the following form 



0. (4.9) 



{-PM + r^) ;! \ +- 1 ^ ^ + 



Qj.2 



dr 



8a2 ypM-r"^ 



-m /(r) = 0. (4.10) 



Setting / (r) = and introducing the tortoise coordinate defined by 



dr 



Eq.(4.10) can be written as a Schrodinger equation 

di^ijj (r) 



dri 



+ [V{r)-n'^] ij{r) = 0, 



where V (r) reads 



V(r) 



M 3r^ PmM mr'^ 

2P^ ^ ~ Al^ ^ 8q;2 ~8a2 



(4.11) 



(4.12) 



Note that the above potential depends only on the separation constant and the Gauss- 
Bonnet coupling constant a. This is the only information that it has from the bulk. It 
docs not depend on the deficit angle /3. This can be understood because b{p) expresses the 
response of the bulk geometry to the presence of the conical singularity, information that 
is entirely encoded in the decoupled equation (4.8), which describes the behaviour of the 
perturbation in the extra-dimensions. We will show in the following subsection that the 
deficit angle explicitly appears in this extra-dimensional equation only if we consider an 



-13- 



angular dependence in the total perturbation as can be seen in Eq.(4.28) together with its 
solution (4.47). 

The above potential is similar to the one found in the calculation of quasinormal modes 
of BTZ black hole [51]. Furthermore, it is exactly the same as the potential we will find 
in the next subsection by solving the Klein-Gordon equation for the background metric, 
provided that we make the identification — )■ The Schrodinger equation (4.11) will 
be solved in the next subsection. The solution implies that 

Si = __|^l + 2^+^l + _j, (4.13) 

where is a positive integer. As this quantity is always negative, it yields a decaying 
solution for the perturbation assuring the stability of the model under scalar metric per- 
turbations. Although in general the mass may also take negative values, the spectrum of 
m can be analysed by looking at the bulk equation (4.8), whose general solution is 

y{p) = Ci 2^1 (^, B; C; q) + 2F,{A': B'; C; q) , (4.14) 

with 



Q = 


cosh^(p/2-\/a) , 


(4.15) 


A = 




(4.16) 


B = 




(4.17) 


C = 


1 - i^4 + 2/2m/a. 


(4.18) 


A' = 


-C/2, 


(4.19) 


B' - 


B-A-C/2, 


(4.20) 


C = 


-2A. 


(4.21) 



For this solution let us choose Ci and C2 such that we can obtain a decaying behaviour which 
ascertains a smooth transition towards infinity. As the arguments of the hypergeometric 
functions need to be real, the first restriction on m comes from the square root: m > 
—2a/l'^. The equality in this equation when substituted in (4.13) corresponds to the strong 
coupling or Chern-Simons limit. Interesting enough, at this limit the extra term disappears 
and we recover the results of a three-dimensional quasinormal modes of the BTZ black 
hole [51]. However, the hypergeometric functions are not well-defined for all the values 
of m above this limit. This strictly depends upon the parameters of each function. A 
careful analysis shows that C2 needs to be zero since it cannot produce decaying solutions. 
Moreover, recalling the properties of these functions, we were able to find a special case 
by setting B = 0, which corresponds to m = 48a /l"^, where the perturbation manifests a 
decaying behaviour. 

We have one more scalar mode namely the hpg mode. As we have already mentioned, 
this mode is coupled with the hte and hro modes. Still using the transverse gauge, it can 
be easily seen, that this scalar mode also decouples from the other two modes and quite 
surprisingly it obeys the same Schrodinger equation as the hoe mode. 
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4.3 Scalar perturbations-solving the Klein-Gordon equation 

In this subsection we will carry out the analysis of the scalar perturbations solving the 
Klein-Gordon equation. We will show that wc get the same results as the ones we got by 
solving the Lichnerowicz equation, but the study of the Klein- Gordon equation gives us a 
better understanding of the stability of the codimension-2 black string. 
Consider the massive Klein-Gordon equation 

= -^dui^f^g^^dN)^ = m'^ . (4.22) 

By using the metric (4.3) and decomposing the scalar field as 

$(i, r, 0, p, e) = Z{t, r, p)S(0)e(^) , (4.23) 
we arrive to the following equation, 

1 d^Z ^ ( ^ 2nh \ dr-Z ^ n^O^ ^ /3f ^^\dpZ_^ d^Z 



p-n? Z \Pr P J Z P Z \ f hj Z Z 

We can apply variable separation method to the 9- and 0-dependent parts of this equation. 
Thus, the solutions to 0(6*) and S(0) are given by 

Q{d) = Ae'^^ + Be-""^ (4.25) 
5(0) = Ce''^ + De-'"'' , (4.26) 

where A, B, C, and D are integration constants while k and e are constants generated 
through the variable separation. The remaining equation can be decoupled in two equations 
with the Ansatz Z{t,r,p) — 'i/{t,r)P{p) as follows, 

dl^ - nM — + 2nn J dr^ - n^dl^ + f + t/^n^ J ^ = (4.27) 

a,^P+(3^ + Qa,P+(-^-m^ + ^)p = 0, (4.28) 

being v'^ the constant of variable separation. Equation (4.27) describes the evolution of the 
scalar perturbation on the brane, while equation (4.28) says how this perturbation behaves 
at different distances away from the brane. 

In order to solve Eq.(4.27) we should rewrite the equation in terms of the tortoise 
coordinate, 

dr 



(4.29) 



Thus, Eq.(4.27) adopts the Schrodinger form. 



0^2 + = V[T{u)]X{t, r,) . (4.30) 
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Here X{t, r*) is a new function defined as 

and the so-called potential term can be written as 

V{r) = — I 2nn - — ) H ^ h v^n^ 



2r 



2rJ 



(4.31) 



(4.32) 



When we compare this potential with the pure 3-dimensional case studied in [51], we see 
that Eq.(4.32) has a u'^n'^ additional correction. This new term connects brane to bulk 
perturbations through the constant u and does not vanish, in general. 

We choose the time dependence as e"*'*'*, so that X{t,r) — e~*'^*it!(r), thus, Eq.(4.30) 
becomes 



+ [u^ -V{r)]R^O. 



(4.33) 



When n(r) is of the BTZ form, n (r) = —M + r /I , the potential (4.32) becomes 



rMe^ 1 — . 



The tortoise coordinate can analytically be found from (4.29) as 

r = -WM coth 



(4.34) 



(4.35) 



Notice that rn — ly M < r < oo corresponds to — oo < r* < 0. 
Let us make the following variable transformation. 



X — 



cosh 



2/rWM 



(4.36) 



with < a; < 1, so that the potential turns into 



V{x) = -\ 

and Eq.(4.33) becomes 

, , (PR , , „ , dR 

4x(l - x)—r + 4 - 6a; — 
dx-^ dx 



4M(1 + i/2;2) -Mx- APe^{x - 1) 
i'^ix - 1) 



0^ 1 + 1/2/2 



X 



X -1 4:{X - 1) 



If 



i? = 0. 



(4.37) 



(4.38) 



In order to write this equation in a more familiar way we make a new substitution as 
follows, 

{x - 1)3/4 



R 



X 



iujl/2VM 



(4.39) 
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Then, Eq.(4.38) turns to be 

—x{x — l)y" + 



iul , 

l-3x + ^=ix-l] 



M 



^ ^ [W ^ 7^ ~ 4M ~ ^ 4(x-l)J ^ " 

(4.40) 

The general solution of this equation can be expressed in terms of hypergeometric 
functions of the second kind 2Fi{a,b; c; x). As we look for stable behaviour, we keep the 
decaying branch of the general solution and establish the restrictions for this stability, 
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y(a;) = (l-x)-(^+^^)/2 2F,ia,b;c;x), 



(4.41) 



where 



b = 



= 1 



1 - Vi + u^P - 
1 - Vi + ly^P - 

iuil 



iuol 
iul 



ih 



s/M 
He 



M 



M 



(4.42) 



Since the (2+1) spacetime is asymptotically AdS, the correct boundary condition to be 
taken in Eq.(4.41) is the flux condition, i.e., 



^ ^ {y*df,y - yd^y*) 



= 0. 



(4.43) 



In order to evaluate this condition at x = 1 we can use the following property of hypergeo- 
metric functions. Given a 2^^1(0', b'\ c'; x), if c' is not a negative integer, the series converges 
when x = 1 if 3fi(c' — a' — b') > 0, and we can write the hypergeometric function as 



2Fi(a',6';c';l) 



T(c')T(c' 



b') 



(4.44) 



r{d - a')r(c' - b') 

The derivative of the hypergeometric functions we are dealing with has the following form, 

ab 



2Fi'{a,b;c;x = 1) = — 2^1(0 + 1, 5 + 1; c + 1; a; = 1) . 



(4.45) 



Using (4.42) we verify that 3?(c- a - 6 - 1) = y/l + u'^P - 1 > since u'^ > 0. Thus, we can 
write (4.45) in terms of F functions as in (4.44) and put back in the flux condition (4.43). 
The new equation is fulfilled when c — a = —N or c — b = —N, where N is a positive 
integer. In this way we obtain the quasinormal frequencies. 



{1 + 2N + Vl + ly^P) . 



(4.46) 



We observe that the imaginary part of these frequencies is negative, thus, displaying the 
stability of the model under scalar perturbations. Notice that in (4.46) a term proportional 
to appears which encodes the information from the bulk and when i/^ — the pure 
three-dimensional BTZ results [51] are recovered. 

To complete our analysis we must solve Eq.(4.28). We have two subcases. 
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1. When /(p) = cosh(/)/2y^) and b{p) — 2^-^sinh(p/2y^), the most general solution 
of Eq.(4.28) is given by 



(4.47) 



/'2z 

where 



a = - I 1 + -r + Vl + 4i/2a - 2 Vl + m^a 



Ci and C2 are constants, and 

d = ^ ^ ~ Vl + 4:iy^a + 2Vl + m'^a^ 

1 / K 

b = - ( 1 + - - Vl + 4z/2q; - 2VI + m^o; 
c = 1 - Vl + 4i/2q; 

Ul + ^ 

2 V /3 

6' = \(^^+'^ + Vl + 4i/2Q; + 2VI + m^o; 
c' = 1 + Vl + 4z/2q; . (4.49) 

A careful study of both hypergeometric functions shows that C2 needs to be zero if we 
want to have a decaying solution. Given a suitable set of parameters the remaining 
hypergeometric function has the pursued behaviour. A special case appears when i/^ 
adopts the following form, 

which produces 6 = 0. This case corresponds to a damped oscillation in p. Therefore, 
the model is well behaved under this kind of perturbation having the quasinormal 
frequencies given in (4.46). 

2. When /(p) = ±1 and 6(p) = 7sinh(p/7), with 7 = {]?■ — ^a)/2, we find the 
following solution for (4.28), 

P{w) = y/2{w - ly''' [Ca 2^1(0, 6; c; w) + C4V2^ 2^^1(0', 6'; c'; «;)] , (4.51) 



where 



Cs and C4 are constants, and 



a = - + - - - Vl + - 4a)(2m2 - 2t/2) 

& - i + I + ^ Vl + (P - 4Q;)(2m2 - 2u^) 
1 

^~ = 2 

a' = i + f + + - 4a)(2m2 - 2i/2) 

c' = ^. (4.53) 

In this case both hypergeometric functions display a growing behaviour. This means 
that the scalar perturbation when going into the bulk propagates without boundary, 
an analogous behaviour to the one found in the five-dimensional black string solution 
of [44]. 

In the last two subsections we studied the scalar perturbations of the five-dimensional 
black string solution of codimcnsion-2. Using both the Lichncrowicz equation and the 
Klein-Gordon equation we found stability under scalar perturbations. For an alert reader 
we comment on the approximation we used. We have chosen the s-wave approximation 
for our ansatz, both for the angular coordinate on the brane and on the transverse space. 
However, on general grounds the Fourier decomposition should also include, the frequency 
modes of the two angular dimensions. Still it is not clear to us whether the deficit angle 
is playing a significant role in the stability analysis. The procedure we have followed in 
the quasi-normal mode analysis shows (see equations (4.27) and (4.28)) that each angular 
coordinate, affects the corresponding spatial dimension. The (f) dependence is reflected in 
the equation for the brane, while the 9 dependence and hence the deficit angle comes with 
the transverse space equation. Any effect of the deficit should be refiected on the solution 
of the transverse piece. However, the only information that equation (4.28) give us is to 
tell us how the perturbation created on the brane is transmitted into the bulk. 

Should we had chosen a frequency mode in the Fourier decomposition for the 9 dimen- 
sion, then we should examine whether the scalar part decouples from the rest of the system 
and even if the variables in the Klein-Gordon equation can be separated. This procedure is 
very involved, due to the presence of the Gauss-Bonnet's quadratic terms and it is beyond 
the scope of the present work. This procedure may give different results, however, even in 
this case, the deficit angle due to its intrinsic geometrical nature, may still has an effect 
only on the transverse space. 

4.4 Vector and tensor perturbations 

We have seen that as far as we are away from the Chern-Simons limit we can always 
apply the transverse and traceless gauge and examine the stability behaviour. We saw that 
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for scalar perturbations the system declines exponentially, meaning that Q is negative, 
even though the mass of the scalar mode takes negative values until a certain limit. This 
feature ensures that our system is stable under scalar metric perturbations. Of course, the 
quadratic nature of the Gauss-Bonnet combination does not give straightforward decoupled 
equations for the scalar part and the decoupling procedure is not easy, as it happens in the 
case of having just the Einstein tensor where the scalar part of the perturbation decouples 
immediately. 

It is known that a category of spherically symmetric vacuum solutions with a Gauss- 
Bonnet term suffer from ghost-like instabihties, and furthermore, there is a strong couphng 
problem close to the Chern-Simons limit [22]. In our case we expect the strong coupling 
problem to be present also in the vector and tensor perturbations. Our analysis has shown 
that under the particular Ansatz that we have chosen for the perturbation there is a pro- 
portionality factor to every equation, namely, (/^ — ia). As ia, all equations are 
identically satisfied. This means that the first order parametrization of the perturbation is 
not significant, and we have to move to the second order. 

We saw at the beginning of this section that the scalar part of our perturbation, can 
be decoupled from the rest of the system, with the help of the transverse and traceless 
gauge. Furthermore, a careful investigation of the perturbation equations unveiled that the 
scalar perturbations do not have any pathological behaviour. To discuss the other kind of 
perturbations we will set these modes equal to zero and examine the rest of the system. 

The equation for the vector mode h^pg is very simple and it is decoupled in a straight 
forward manner from the rest of the system 



The solution of this equation is trivial and it gives information only about the transverse 

space. As we can see there is no information about the branc. Before making a comment 
on that let us sec the equation for hte and hro modes. Making use of the transverse gauge, 
these modes decouple from each other and the result is that both of them obey the same 
equation as h^po, namely equation (4.54). This is quite unnatural from many aspects. First 
of all since the brane dependence is arbitrary, this means that it could be satisfied by a 
mode that has positive fl, signaling this way departure from stability. Furthermore, the 
system is degenerate since these three modes behave exactly the same, while at the same 
time we can not make any solid prediction about its brane behaviour. This could be a 
signal of a strong coupling problem, but it is not clear whether this is exactly the case or 
not. 

We still have the tensor part plus one more vector. It can be seen that choosing any of 
the remaining perturbation equations, with the help of the transverse and traceless gauge, 
we can set them identically equal to zero. This is again an unnatural result. Either it 
means that we have no propagating degrees of freedom, which indeed could be the case 
because in three dimensions we do not have a graviton, or again we are in a strong coupling 
regime, which restricts us from making any conclusion from the linear terms and we have 
to move forward to higher order ones ^. However, the case where we have no propagating 

^In [52] the holographic description of the Gauss-Bonnet theory in five dimensions was discussed and it 




(4.54) 
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degrees of freedom seems more probable, since in the s-wave approximation that we have 
chosen, the brane can "see" the transverse space through the angular symmetry that is 
manifested both on the brane and on the transverse space. The absence of a graviton mode 
on the brane, due to its three dimensional nature is reflected also in transverse space due 
to the aforementioned symmetry. 

5 Conclusions 

We studied the metric perturbations of a five-dimensional black string in codimension-2 
braneworlds with a Gauss-Bonnet term in the bulk. After reviewing the general formal- 
ism of linear metric perturbations, we derived the modified Lichnerowicz equation in the 
presence of the Gauss-Bonnet term. As an application, we considered the six-dimensional 
spherically symmetric Gauss-Bonnet black hole solutions and we applied the modified Lich- 
nerowicz equation to these solutions to derive the known results for the tensor perturba- 
tions. 

We considered the scalar perturbations of a five-dimensional black string solution of 
codimension-2 braneworlds with a Gauss-Bonnet term. We carried out the analysis using 
both the Lichnerowicz equation and solving explicitly the Klein-Gordon equation. Away 
from the Chern-Simons limit, we showed that the results from the two methods coincide. 
The behaviour of the black string under scalar perturbations can be described by two 
equations. One of them describes the evolution of the scalar perturbation on the brane 
while the other equation shows how this perturbation behaves at different distances from 
the brane. The evolution of the black string on the brane is controlled by the quasinormal 
modes which are similar to the quasinormal modes of the BTZ black hole with the addition 
of an extra term which has the information from the bulk. The stability analysis shows 
that the black string is stable under scalar perturbations. 

We also studied the vector and tensor perturbations of the black string solution. We 
found that for the vector modes we can only get information on their behaviour in the bulk 
while for the tensor modes we did not find any physical propagating modes. We mainly 
attribute this behaviour to the specific symmetries of the considered black string solution 
but also may be an indication of a strong coupling problem signaling the need to go to the 
next order in perturbation theory. 

It is interesting to extent the analysis to the six-dimensional black string solutions 
derived in [40]. In six dimensions there is also a limit where the theory becomes strongly 
coupled. Away from this limit, we expect the application of the modified Lichnerowicz 
equation to this solution to simplify significantly the calculations but the main problem 
which has to be addressed is the presence of matter in the bulk. The matter is necessary 
to support the black string solution on the brane and its extension into the bulk. However, 
there is a limit where the matter in the bulk decouples and the general formalism developed 
in this work can be applied. This issue is under investigation. 



was found that there is a particular Weyl anomaly that prevents the Gauss-Bonnet theory to go smoothly 
to the Chern-Simons limit signaling the breakdown of perturbation theory at linear order. 
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A The Lichnerowicz equation for the scalar modes 

In this appendix we give some technical details for calculating the scalar modes. First, 
we take the traceless equation, solve it for the function hrr, and then substitute into the (99) 
component of the Lichnerowicz equation (2.32). Also wc substitute the function h,.,- into 
V^h^p equation, from which we solve for dj-hrp- Then, we replace both drhrp and drdphrp 
into the (66) equation. We follow the same steps also for the (pp) equation of (2.32). From 



tanh^ ( Yy— i , ^ 

the resulting two equations we take the combination Aa/S ) \ [PP) — (^^) and we find 

the following equation for the two scalar modes 



(2 Pan'- {PM - r') cosh (^) sech^ (^) ) 

4/2 {pM-r'^) a 

^ (3r2 - P (3M - a^f)) (l - % cosh (^^) ) - 5 [PM - r') cosh 

7 ^ hpp 
8P {P M - r2) sech"^ (^j 

6csch-^ (^) sinh^ (^) ^/32 (-3 + 4 cosh (^)) sinh (^) 

7-) — -F= 8p hee 7 N (Jphpp 

^^cosh^l^fe) 

dpphee 2al3Hanh'' (^^^ dpphpp 



P (l + cosh(^)' 

(PM 3r2) sech^ d,heo O^M r') sech^ clrhe 

32 {PM - 3r2) aP^sch"" (^) sinh*^ (^7^) ^rhpp 



32 {PM - r2) a^^sch' (^) sinh^ (^^j 9,, V 

— (A.l) 



/4 



We need one more equation in order to find a solution for the scalar perturbations. The 
second equation comes from the {tp) component of (2.32), and it can be derived in a similar 
fashion as (A.l). Again we solve the traceless condition for and substitute into the {tp) 
equation. We also replace this equation into V^hAp and solve for drhrp- We do the same 
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1 + 7cosh (^)] scch^ (4^) tanh^ (^^ 



Qpp'^ {r, p) 



dpppU (r, p) 



2 (Z^M - 3r^) asech^ (^g) tanh^ (^g 
2(m-3r^)asech^ (2^) (275) 



drU (r, p) 



/4 



9^^!* (r, p) 



128 (/2M - 3r2) at csch^ (^^) sinh^ 



128 (/2M - 3r2) csch^ (^^) sinh^ (2^) 

Z4 



(r, p) 



drrpU (r, p) = 



(A.5) 



After some algebra the above equation takes the following elegant form, 

3 



d 



9p \ {PM - r2) cosh (^275) 



u (r, p) + 



1 - 2 cosh^ ( ^ 

2 V" 



-2 cosh 



2# 



p \ (9^ M (r, p) ^9 



r fr^M — r 



du{r,p) 



2^ J ap2 



r cosh 



= 



(A.6) 



This equation is fully integrable in the p coordinate. By setting the integration constant to 
zero and separating the variables with the Ansatz u (r, p) — f (r) y (p) , we get the following 
two differential equations, 



with m being the variable separation constant. 



p \ d'^y{p) 



\foL ) d p^ 



+ my{p)^0, 
(A.7) 



0. 
(A. 
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